Anomalous Exciton-Condensation in Graphene Bilayers 
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In ordinary semiconductor bilayers, exciton condensates appear at total Landau level filling factor 
vt = L We predict that similar states will occur in Bernal stacked graphene bilayers at many 
non-zero integer filling factors. For vt = —3, 1 we find that the superfluid density of the exciton 
condensate vanishes and that a finite-temperature fluctuation induced first order isotropic-smectic 
phase transition occurs when the layer densities are not balanced. These anomalous properties of 
bilayer graphene exciton condensates are due to the degeneracy of Landau levels with n = and 
n — 1 orbital character. 
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Introduction — When two semiconductor quantum wells 
occupied by half-filled Landau levels are narrowly sep- 
arated, the bilayer system ground state spontaneously 
establishes inter-layer coherence These broken sym- 
metry states possess a condensate of pairs (each com- 
posed of an electron in one layer and a hole in the other) 
which opens up an energy gap responsible for a total 
filling factor vt — 1 quantum Hall effect, supports a dis- 
sipationless counterflow (excitonic) supercurrent, and is 
responsible for a wide variety of incompletely understood 
transport anomalies 0,13. Recently an interesting new 
type of bilayer two-dimensional electron system has be- 
come available (3, [3 which consists of two carbon-atom 
honeycomb-lattice (graphene) layers separated by a frac- 
tion of a nanometer. The electronic structure of graphene 
bilayers is full of surprises [3 because of an interplay be- 
tween the sublattice pseudospin chirality of each layer Q 
and the Bernal stacking arrangement, particularly so in 
the presence of an external magnetic field. In this struc- 
ture, one of the two-carbon atom sites in each layer has a 
near-neighbor in the other layer and one does not. Inter- 
layer hopping drives electrons on the closely spaced sites 
away from the Fermi level, leaving one low-energy site 
for carbon 7r-orbitals in each layer. Because the hopping 
process between low-energy sites occurs in two steps, via 
an intermediate high-energy state, it turns out[6| that in 
the presence of a magnetic field both n = and n = 1 
Landau level orbitals have zero kinetic energy, and that 
the corresponding wavefunctions are localized in oppo- 
site layers for opposite (K and K') graphene valleys [5(- 
Even though the two-layers are close together and band 
eigenstates at zero-field are a coherent combinations of 
individual layer components, the zero-energy strong field 
states occupy a definite layer. 

In the absence of electron-electron and Zeeman in- 
teractions, a group of eight degenerate Landau levels 
(the bilayer octet) develops at the Fermi level of a neu- 
tral graphene bilayer because of the degeneracy between 
n = and n = 1 states combined with spin and val- 
ley degeneracy. The octet is revealed by a jump in the 



quantized Hall conductivity [3 from —A(e 2 /h) to 4(e 2 //i) 
when the charge density is tuned across neutrality in 
moderately disordered samples. When the sample qual- 
ity is sufficient for interactions to dominate over disorder, 
quantum Hall effects occur [H, 0, [13 at all intermedi- 
ate integers. Spontaneous valley-symmetry breaking is 
expected[9( in all but the Vt = (half-filled octet) case. 
Interactions favor states with spontaneous coherence be- 
tween valleys, and hence between layers, because this 
broken symmetry does not require charge transfer be- 
tween layers. In this Letter we illustrate the unique and 
rich properties of bilayer graphene exciton condensates 
by concentrating on the simplest case in which the octet 
is occupied either by single majority-spin (i>t = —3) or a 
single minority-spin {vt — 1) Landau level, and allowing 
for an external potential difference Ay between the lay- 
ers. We demonstrate that the superfluid density vanishes 
at these filling factors, and that a finite-temperature fluc- 
tuation induced first order isotropic-smectic phase tran- 
sition occurs when the layer densities are not balanced. 
These anomalous properties of bilayer graphene exciton 
condensates are due to the degeneracy of Landau levels 
with n = and n — 1 orbital character. Below we first 
discuss the mean-field ground state, demonstrating that 
it exhibits spontaneous inter-layer phase coherence for 
small Ay at the filling factors of interest. We then ex- 
plain why the superfluid density vanishes and discuss two 
important consequences, namely that the phonon collec 



tive mode has quadratic rather than the expected 11 1 
linear dispersion and that long-wavelength instabilities 
appear when Ay =/= 0. We conclude with some specula- 
tions on the experimental implications of our findings. 

Octet mean-field theory for unbalanced bilayers — Bilayer 
graphene's N = Landau-level octet is the direct prod- 
uct of three S — 1/2 doublets: real spin and which- 
layer pseudospin as in a semiconductor bilayer, and the 
Landau-level (n — 0,1) pseudospin which is responsible 
for the new physics we discuss in this paper. The octet 
degeneracy is lifted by the Zeeman coupling A z , assumed 
here to maximize spin-polarization, and by Ay which 
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is denned so that it favors top layer occupation when 
positive. In graphene bilayers, Ay also drives @, Q a 
small splitting of the Landau-level pseudospin that plays 
a key role in this paper: — aAyfiw c /7i where 

hw c = 2.14£?[Tesla] meV, 71 ~ 400 meV is the inter- 
layer hopping energy, and a = +1(— 1) for top (bottom) 
layers. The fact that LL 1 has a higher energy than LL 
in the top layer and a lower energy in the bottom layer 
will prove to be important. 

The octet mean-field Hamiltonian can be separated 
into single-particle and exchange contributions: 
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(repeated indices are summed over) where the single par- 
ticle energy (which includes the Hartree capacitive term) 
is 

lAz Ay d (v \ 

E a ,s,n = -s— a— + anA LL + — y - - h^j , (2) 



and the exchange interactions 
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In these equations, <? is the Landau level degeneracy, 
d = 3.337A is the interlayer separation, l B is the mag- 
netic length, Az = gp. B B is the Zeeman coupling, 
n = 0, 1 are LL indices, a, (3 = t,b for top(bottom) 
layers, r, s = 1(— 1) for up(down) spins and t — b,b — 
t. All energies are in units of e 2 /el B . The average 
value of the operators p^nv,t,m = Ex c l, s ,ni,x c 0,t,n 2 ,x 
(where s n x creates an electron with guiding-center 
X in the Landau gauge and layer a, spin s and Lan- 
dau level character n) must be determined self consis- 
tently by occupying the lowest energy eigenvectors of 
H HF . The form factors (F 00 (q) = e^ 1 ^ / 4 , Fi (q) = 
{iq x + q y )l B e- {qlB)2/ yV2 = [-fbl(-q)]* and F n (q) = 
(1 - (ql B ) 2 /2)e^ qlB ^^) reflect the character of the two 
different quantum cyclotron orbits. 

Although an infinitesimal Ay would be sufficient to 
produce complete layer polarization in a non-interacting 
system, a finite value is required once interactions are 
accounted for. Spontaneous interlayer phase coherence 
arises because it is able to lower energy by inducing a gap 
at the Fermi level even when both layers are partially oc- 
cupied. When the layer index is viewed as a pseudospin 
[H, the phase-coherent state corresponds to an x — y 
easy-plane ferromagnet and the interlayer phase differ- 
ence corresponds to the azimuthal magnetization orien- 
tation. In this language Ay is a hard-direction external 



field which gradually tilts the magnetization direction to- 
ward the z direction. We find that for Ay > Ay (B) , 
with A* v (B) /(e 2 /el B ) ~ 0.001 at B = 10T, the system 
exhibits full charge imbalance with all electrons (all mi- 
nority spin electrons in the Vt — 1 case) occupying one 
layer. The critical Ay is given by 
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(4) 

In the graphene bilayer case, the Ay sufficient to achieve 
full pseudospin polarization is reduced compared to the 
semiconductor case because the layers are close together. 
For Ay < Ay, the charge-unbalanced mean-field ground 
state consists of a full Landau level of states [t| [13, EH 
which share partially polarized layer and n — Landau- 
level pseudospinors: 
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10), 
(5) 

where we have dropped the irrelevant spin degree-of- 
freedom. In Eq. §5§ , tp is the spontaneous interlayer phase 
and cos (0v) = Ay / A* v (B) . 

Vanishing superfluid density at zero bias — The energy 
cost of phase gradients is the key property of any su- 
perfluid. In normal superfluids, including semiconductor 
bilayer exciton superfluids, the leading term in a phase- 
gradient expansion of the energy-density is proportional 
to |V</?| 2 . We now show that the coefficient which spec- 
ifies the size of this term, often referred to as the super- 
fluid density, is zero in bilayer graphene exciton conden- 
sates at zero bias. We proceed by explicitly constructing 
a pseudospin wave state in which Vip = qx is constant: 
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where A q x — u q c\ x + v q c\ 1 v and u;. 



(6) 

\,x - u q^b,o,x t u Q u b,i,x allu d q + v q = 1 - In 
this wavefunction the factor exp(iqX) is responsible for 

the phase gradient. The term proportional to v q in Eq. 
([6|) allows for the crucial possibility, unique to bilayer 
graphene, of combining the phase gradient with an ad- 
mixture of the n = 1 wavefunction. 

In mean-field theory the energy is proportional 
to the square of the density-matrix. We now show 
that it is possible to choose v q so that the density 
matrix is unchanged to first order in q and that the 
superfluid density vanishes as a consequence. Since 
v q must vanish for q — > to reproduce l^o); we have 
v q rsj q and u q ~ 1 to this order. It is then easy to 
see that the density-matrix component within each 
layer is unchanged to leading order in q. In quantum 
Hall bilayer exciton condensates, the superfluid den- 
sity is due to reduced interlayer exchange energy in 
the presence of a phase gradient 0, and hence to 
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FIG. 1: (Color online) Collective mode dispersion u q for bal- 
anced Octet QHF in units of interaction strength e 2 1 els = 
11.2(B[T]) 1/2 meV as a function of qls at a magnetic field 
of B = 10 Tesla. The green (circle) curve shows the in- 
version symmetric n = to n = 1 transition mode with 
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dispersion which is unrelated to superfiuid 



properties. The red (box) curve shows the inversion asym- 
metric superfiuid Goldstone mode with the anomalous long- 
wavelength dispersion lim 9 ^o ~ q 2 . The blue (triangle) 
curve is a gapped mode not essential for our discussion (see 
text for details). 



changes in the interlayer density matrix, p t b(r,r') 
*J(r)*b(r')\ where the field operator *t 6) (r) 
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with 4> n {x) the harmonic oscillator state with 
orbital Landau character n. At zero bias, 

Eq. ||SJ) with tp = gives pf b s (r, r') = 
l/(2i») Ex i x - x ) ^o (as' - X) e - iX (y-v')/el wh ii e 
Eq. (JBJ) implies 

pf fc ^(r,r') = ^^^ (a; -X) e -^(^')A 2 B e ^ 
v x 

x [u g o (x' - X) + v q <h. {x' - X)] (7) 

Since <fii(x) = ^/2{x/Ib)4'o{x) and q(x' — X) is small be- 
cause of the localized oscillator wavefunctions, it follows 
that pf b s (r, r') is altered only by an irrelevant phase fac- 



tor to first order in q if we choose u„ 



The leading change in \ptb\ is therefore proportional to 
q 2 and the leading energy change proportional to q*. 
Unusual collective mode dispersion — We now consider 
collective excitations of the bilayer quantum Hall exciton 
condensate at vt = —3, 1, first for balanced (Ay = 0) 
bilayers. As discussed above, the ground state is a full 
Landau level with shared layer symmetric n = pseu- 
dospin orbitals | + ,0) = 0) + |6, 0)) /a/2. The pseu- 
dospinor fluctuations that are quantized in the system's 
collective modes mix in components from the three or- 
thogonal pseudospin states: | + , 1) also symmetric in layer 



indices and | — , 0) and | — , 1) antisymmetric in layer in- 
dices. Fig. [1] shows the collective mode spectrum calcu- 
lated in the time-dependent mean-field theory [la ]. The 
inversion symmetric mode (transition | + ,0) — > |+, 1)) is 
unrelated to superfluidity and has been discussed previ- 
ously Q. For q = the two asymmetric modes corre- 
spond respectively to global rotations of the |+, 0) pseu- 
dospinor toward the states | — ,0) and | — , 1) respectively. 
The later rotation costs a finite energy because the ex- 
change energy in n = 1 is smaller than in n = 0. 

One physically transparent way of performing these 
collective mode calculations is to construct a fluctuation 
action in which each transition has canonically conjugate 
density p and phase ip components corresponding to the 
real and imaginary parts of the final state component in 
the fluctuating spinor. The fluctuation action 



(8) 



S[p, ip] = S B [p, (p}~ J duj J d 2 q £[p, ip], 
contains a Berry phase term Sb flfil ] 
S B [p,(p] = J duj J d 2 q ipL q XV>q - <pL q Z>Vq 



(9) 



where p q = (pi, q , p 2 ,q, P3,q) and V = -M 3X 3, and an 
energy functional £ [p, <p] closely related to the discussion 
in the preceding section 



(10) 



Here, A p (q) and A ip (q) capture the energy cost of small 
pseudospinor fluctuations and can be evaluated explic- 
itly. Because p fluctuations change the charge distribu- 
tion in the system they remain finite for q — > and do 
not play an essential role in our discussion. At long wave- 
lengths, we find that the inversion asymmetric block in 
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with "• • •" representing terms higher order in qls- 

In a semiconductor bilayer only n — phase fluc- 
tuations are possible without paying a kinetic energy 
penalty. In the strong magnetic field limit the super- 
fluid density is therefore proportional to the interlayer 
exchange-interaction constant -^iqoi?- For the bilayer 
graphene octet, however, the energy cost of phase vari- 
ation is proportional to the smallest eigenvalue of the 
matrix in Eq. (fTTj) and this has the g 4 -long wavelength 
dependence anticipated above. Indeed the eigenvector of 
this quadratic form captures the orbital character of the 
Goldstone mode 
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corresponding to the v q /u q ratio in that analysis. We 
find that the leading small q behavior for the energy func- 
tional is 



£_[V>_]=/3(?J B )V.+ 



(13) 



where j3 = 0.093(e 2 /els) at 10T and ip_ is the ampli- 
tude of the lowest-energy eigenmode of the phase matrix 
A v (q). Because of the conjugate relationship between p 
and <p fluctuations, the collective mode energy is propor- 
tional to the square root of the eigenvalues of A v and A p 
so that the quadratic Goldstone mode dispersion simply 
signals the vanishing superfluid density. 
Long-wavelength instability at finite bias — The q A in- 
teraction energy cost of inter-layer phase gradients holds 
for balanced and unbalanced bilayers. In unbalanced lay- 
ers, however, the Landau level splitting All means that 
the single-particle energy is lowered by transitions to the 
n = 1 Landau level in the bottom layer. For unbalanced 
layers we find that the phase fluctuation action goes to 
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The energy cost of phase fluctuations is again propor- 
tional to the smallest eigenvalue of the phase matrix, 
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which is negative at small q, indicating that a uniform 
condensate is unstable when All 0. The energy func- 
tional in Eq. ([T5|) reduces to the classical Swift-Hohenberg 
(SH) model [17] Hamiltonian: 
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when we set A = -(cos (6 V ) A LL ) 2 /8/3, $ = VWl 2 B 
and (q^ls) 2 — cos {By) All/^P- Using the detailed mi- 
croscopic form of the mean-field energy functional, we 
estimate that A ~ (<7o's)( e2 / e ^)- The SH model ex- 
hibits a fluctuation-induced first order smectic-isotropic 
phase transition from a stripe ordered phase with (</?-) = 
Acos(qo -r) to a disordered phase with (<p~) — 0. Follow- 
ing the self-consistent Hartree approximation analysis [l7| 
of the Swift-Hohenberg model, for Ay < Ay we estimate 
the transition temperature to be 
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which is typically below lOmK. 

Discussion — We identify the T > T c phase with a nor- 
mal quantum Hall ferromagnet and the T < T c phase 
with a quantum Hall smectic state which should exhibit 
anisotropic transport properties. For T > T c we expect 



properties similar to those observed in semiconductor 
bilayers [2| except that the coherent interlayer tunneling 
processes, which plays a key role in tunneling experi- 
ments, should be essentially absent. In graphene bilayers 
therefore, spontaneous coherence is likely to be most con- 
veniently manifested by strongly enhanced Hall drag. Fi- 
nally we note that trigonal warping terms, which we have 
neglected, will reduce A£l[9|] by at most 5% for magnetic 
fields of interest, that correlation effects we have not con- 
sidered could contribute positively or negatively to the 
superfluid density, and that the small superfluid densities 
in this system might lead to important thermal fluctua- 
tion effects beyond those considered here. This work was 
supported by the NSF under grant DMR-0606489 (AHM) 
a grant by NSERC (RC) and the State of Florida (YB). 
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